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We develop a simple model for the evolution of a nucleon spin in a hot and dense nuclear medium. 

A given nucleon is limited to one-dimensional motion in a distribution of external, spin-dependent 
scattering potentials. We calculate the nucleon spin autocorrelation function numerically for a 
variety of potential densities and distributions which are meant to bracket realistic conditions in a 
supernova core. For all plausible configurations the width T of the spin-density structure function is 
found to be less than the temperature T. This is in contrast with a naive perturbative calculation 
based on the one-pion exchange potential which overestimates T and thus suggests a large suppression 
of the neutrino opacities by nucleon spin fluctuations. Our results suggest that it may be justified to 
neglect the collisional broadening of the spin-density structure function for the purpose of estimating 
the neutrino opacities in the deep inner core of a supernova. On the other hand, we find no indication 
that processes such as axion or neutrino pair emission, which depend on nucleon spin fluctuations, 
are substantially suppressed beyond the multiple-scattering effect already discussed in the literature. 

Aside from these practical conclusions, our model reveals a number of interesting and unexpected 
insights. For example, the spin-relaxation rate saturates with increasing potential strength only if 
bound states are not allowed to form by including a repulsive core. There is no saturation with 
increasing density of scattering potentials until localized eigenstates of energy begin to form. 

PACS numbers: 26.50.+x, 13.15.+g 


I. INTRODUCTION 

The numerical modeling of stellar collapse and super¬ 
nova explosions has enjoyed enormous progress which is 
not quite matched by the poor development of the mi¬ 
croscopic input physics. This unsatisfactory state of af¬ 
fairs has recently motivated a number of authors @1 
to turn their efforts toward a better understanding of 
one particular key ingredient, the neutrino opacities in a 
hot and dense nuclear medium, for which only crude pre¬ 
scriptions have thus far been used in supernova codes. 
The emphasis of these works has been to address the 
composition dependence of the opacities (hyperons may 
be important), and to include many-body effects which 
range from Pauli phase-space blocking over the inclusion 
of modified nucleon dispersion relations to “screening,” 
i.e. nucleon-nucleon correlations. 

We study a further aspect of in-medium neutrino re¬ 
action rates which was not mentioned in these papers, 
the topic of nucleon multiple-scattering effects. Their 
possible importance derives from the observation that 
the scattering of neutrinos on nonrelativistic nucleons is 
dominated by axial-vector current interactions, i.e. it is 
primarily a spin-dependent phenomenon. The nucleon- 
nucleon interaction contains a strong spin-dependent 
piece, leading to a large nucleon spin-fluctuation rate in 
a hot nuclear medium. In the nonrelativistic limit, these 
spin fluctuations are the sole microscopic source for the 


emission of neutrino pairs or axions from a supernova 
core or a neutron star |Q||, processes which are pertur- 
batively described as nucleonic bremsstrahlung. More¬ 
over, spin fluctuations inevitably lead to a reduction of 
the neutrino-nucleon interaction rate and to an in¬ 
crease of the energy transfer rate between the neutrino 
and nucleon fluids ||||. 

The key quantity for our discussion is the rate T of 
spontaneous nucleon spin fluctuations, or equivalently, 
the rate by which a nuclear spin in the medium loses 
memory of its orientation. Of course, this process is de¬ 
scribed by a single constant T only if the polarization of 
a given spin decays exponentially as e~ rt ; in general the 
temporal evolution of the spin autocorrelation is more 
complicated. 

A perturbative calculation of, say, neutrino pair emis¬ 
sion as a bremsstrahlung process NN —> NNvv by virtue 
of a one-pion exchange potential implies a value for T 
much larger than a typical temperature T of a supernova 
core. On the other hand, the resulting suppression of 
the neutrino opacities is in blunt conflict with the large 
duration of the SN 1987A neutrino signal unless one dis¬ 
cards as background all of the late-time events |J. It 
was thus concluded that the true T in a real SN core is 
probably not much larger than T. Additional theoretical 
support for this empirical conclusion arose from a mod¬ 
ified /-sum rule for the spin-density structure function 
which was derived by one of us [^0| . 
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We presently study a one-dimensional toy model for 
the nucleon spin relaxation in a “background medium” 
of spin-dependent potentials. We believe that our model 
captures the essential physics of the nucleon spin evo¬ 
lution, yet is simple enough to allow for an intuitive un¬ 
derstanding and a numerical non-perturbative treatment. 
The conclusion will be that for plausible configurations 
of density and potential strengths we always have T < T, 
suggesting that multiple-scattering effects are indeed not 
a dominating feature of the neutrino opacities in the deep 
inner core of a supernova. Of course, it also means that 
the rates for axion or neutrino pair emission are much 
smaller than given by a naive bremsstrahlung treatment. 
Ironically, then, precisely because perturbation theory is 
bad in a SN core implies that a naive calculation of vN 
scattering, ignoring nucleon spin fluctuations, is a much 
better approximation than it deserves! 

In Sec. II we describe and justify the physical structure 
of our model. In Sec. Ill we discuss our technique for 
the numerical calculation of the spin-density structure 
function. Sec. IV is devoted to numerical results, and in 
Sec. V we conclude with a discussion and summary. We 
use natural units, h = c = ks = 1, throughout. 


II. MODEL 

A. Spin-Density Structure Function 

As a starting point for our model construction we con¬ 
sider a nonrelativistic single-species medium consisting of 
nucleons N. In this case the main quantity of interest for 
neutrino opacities and axion emissivities is the positive 
definite spin-density structure function 

/ +oo 

dt e iuJt R(t, k). (1) 

- OO 

It is the Fourier transform of the normalized spin auto¬ 
correlation function 

R(t,k) = -^(ff(f,k) • <r(0,-k)), (2) 

OTijV 

where <t(cc) = ( x ) <T ' t / ; ( x ) ' s the spin-density opera¬ 

tor with ip(x) the nonrelativistic nucleon field operator, 
a Pauli two-spinor, while er is the vector of Pauli matri¬ 
ces. Further, njy is the nucleon number density and (■ ■ •) 
denotes a thermal ensemble average. The normalization 
factor is given by (<r(0,k) • <r(0, — k)) = s(s + 1) njv for 
all k with s = ^ the length of the spin. Therefore, if the 
nucleon spins are uncorrelated with each other we have 
l?(0,k) = 1, corresponding to 

/ + OO 7 

.^ s <“- k > = i- < 3 > 

We are primarily interested in the relaxation behavior of 
R(t, k). Therefore, we may equally imagine that possi¬ 
ble correlation effects are absorbed in the definition of R 


so that f?(0, k) = 1 can be imposed as a normalization 
condition. 

While S(u>, k) is what appears as a scattering kernel in 
the Boltzmann collision equation for neutrinos, it is not 
the most practical quantity for a theoretical discussion. 
It obeys the principle of detailed balancing, 

S(-u,k) = e~ u / T S(u,K), (4) 

and thus is not symmetric. Correspondingly, R(t, k) has 
a nontrivial imaginary part. However, it is enough to 
consider a symmetric correlator 

1 ,'\ _ 3 (<x(i, k) • <x(0, -k) + <x(0, k) • &(t, -k)} 

k) — 


which is real because R(t, k) = \[R(t, k) + R*(t, k)] = 
k) +R(—t, k)] and thus is the closest quantum ana¬ 
logue to a classical correlation function. It is equivalent 
to a structure function which is even in u>, 

( 6 ) 


The normalization remains unchanged. We can always 
recover the scattering kernel by virtue of 


S(u), k) 


25V, k) 

1 + e~ ul / T ’ 


(7) 


bringing out detailed balance explicitly. 


B. Long-Wavelength Limit 

In the absence of spin-dependent forces the evolution of 
the nucleon spins is trivial. In this case the w-dependence 
of S(uj, k) is exclusively caused by nucleon recoils, im¬ 
plying 5(w,0) = 2ttS(uj). For a nonvanishing momentum 
transfer k, the width of 5 is given by the nucleon recoil 
and thus is suppressed by the inverse nucleon mass. 

The presence of spin-spin interactions changes this pic¬ 
ture in that there will be correlation effects. If we ignore 
the tensor force as in Ref. (Q, the total nucleon spin is 
still conserved in binary collisions, allowing for spin dis¬ 
sipation only by diffusion. While the structure of S(u), k) 
is now more complicated, it still has nonvanishing power 
only for w 2 — k 2 < 0. Put another way, the “screening 
corrections” of Ref. || primarily lead to a reshuffling of 
the power of S(w, k), but it continues to have power only 
where it has already from recoil effects. 

We are primarily interested in the opposite situation, 
i.e. the effect of the tensor force which allows for local 
spin dissipation by coupling the nucleon spins to the or¬ 
bital motion. It creates power in S(u>, k) for w 2 — k 2 > 0 
and, in particular, has a nonvanishing width even in the 
long-wavelength limit k —> 0. Among other consequences 


2 







this implies that neutrino pairs and axions can be emit¬ 
ted or absorbed, processes which are perturbatively de¬ 
scribed as bremsstrahlung NN «-> NNvv. Without the 
tensor force such processes do not occur as can be shown, 
for example, by an explicit calculation. The bremsstrah¬ 
lung picture suggests that S(oj, k) depends only weakly 
on k, essentially because momentum conservation can be 
achieved among the nucleons alone. Therefore, it is prob¬ 
ably enough to develop an understanding of the behavior 
of the structure function in the long-wavelength limit, 
S(u)) = linik-^o S(u>, k). 

Our simple final approach is, therefore, to study the 
single-nucleon spin autocorrelation function as a proxy 
for the collisional broadening of the full dynamical spin- 
density structure function. 

If the effect of collisions on the spin evolution of a given 
nucleon can be pictured as a sequence of uncorrelated 
“kicks,” then the autocorrelation function must decay 
exponentially, 

£(i)=e- r|t l, (8) 

where the quantity T is what we call the spin-fluctuation 
or spin-relaxation rate. If the spin is completely random¬ 
ized in a given collision, one can show that T is identical 
with the collision rate T co n. The structure function cor¬ 
responding to Eq. (||) is a Lorentzian, 

2 F 

SM = - T -- s . (9) 

A (classical) bremsstrahlung calculation involving single 
“kicks” yields S(u) = 2 T/u 2 . The Lorentzian structure, 
which we obtain directly from a Fourier transform of 
Eq. (||), is equivalent to a resummation of a sequence 
of uncorrelated kicks, i.e. it arises from the multiple¬ 
scattering nature of the spin motion caused by a random 
sequence of collisions. 


C. Estimate of T 

The perturbative estimate of T is based on a cal¬ 
culation of the bremsstrahlung process NN —> NNvv 
(Appendix A). Using a one-pion exchange potential and 
nondegenerate nucleons it is found to be approximately 

T/T = 1.24/514 T 3 o 1/2 where pu = p/10 14 g cm -3 and 
T 30 = T/30 MeV. T/T of a few is enough to imply a sig¬ 
nificant neutrino cross-section suppression (Appendix B). 

The question if the perturbatively estimated T is re¬ 
alistic can be addressed by a simple classical estimate. 
The nucleon spin relaxation is caused by collisions with 
other nucleons. The most effective relaxation is achieved 
if the spin is completely randomized in each collision so 
that T = T co n. The maximum collision rate, on the other 
hand, is estimated by the time it takes a given nucleon 
to cross the average distance to the next scattering cen¬ 
ter. For nondegenerate nucleons in thermal equilibrium 


we have ^m^v 2 = § T, leading to a typical velocity of 
v = (3T/toat) 1//2 = 0.31 TpJ 2 . A typical inverse distance 
between nucleons is n]{ 3 = 77.3 p\{ 3 MeV. At this den¬ 
sity the pion-fields produced by the nucleons overlap so 
that the average distance between nucleons should, in¬ 
deed, represent the distance between scattering events. 
Therefore, we estimate that the spin-relaxation rate can¬ 
not be larger than 

T-l 1/3 

» -^r- = 0.8 p\i 3 T 3 - 0 1/2 , (10) 

an estimate which we think is a generous upper limit be¬ 
cause a full spin randomization in each collision is by no 
means assured. Therefore, the perturbative result looks 
too large even for a density 10 14 g cm -3 , not to mention 
supranuclear densities, and the spin relaxation rate is ex¬ 
pected to be significantly smaller than the temperature. 

D. One-Dimensional Model 

The way we performed the classical estimate of the re¬ 
laxation rate essentially defines our model. We picture 
the spin to evolve due to spin-dependent potentials which 
the nucleon encounters as it moves. The spin evolution 
is thus determined by the potential seen along the nu¬ 
cleon path which we picture as being one-dimensional— 
the nucleon motion in coordinate space should not have a 
decisive impact on the (three-dimensional) spin motion. 

We are thus led to ask the following question. If we 
consider a nucleon, constrained to one-dimensional mo¬ 
tion, how does the relaxation behavior of the spin depend 
on the density and strength of the scattering potentials? 
Can the relaxation rate become as large or even larger 
than the temperature for configurations which mimic the 
situation in a supernova core? 

To address these questions we calculate the spin cor¬ 
relation function for a nucleon which is confined to a 
one-dimensional box of length L. We place a number of 
scattering potentials into this box at random locations, 
calculate numerically the exact R(t) or S(u>), and then 
take an ensemble average over many configurations. 

A scattering potential at location xo is taken to in¬ 
teract with the test nucleon according to a Hamiltonian 
density of the form 

Hint = t*{x - Xq) ■ <t{x). (11) 

The classical “magnetic dipole” is 

H(x) = m V(x) (12) 

with m a unit vector pointing in some random direction 
while V (x) is a function with the dimension energy. We 
model it by a Gaussian with 

V(x) = V 0 e - ( x/6)2 (13) 
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where b is a length scale. 

In order to avoid bound states we will sometimes add a 
repulsive diagonal potential so that the total interaction 
becomes 

=V(x- x 0 ) ^ip\x )(m ■ a + l)ip(x). (14) 

For nucleons polarized in the m-direction this implies 
that a spin-down neutron feels no potential at all while 
a spin-up neutron feels twice the (repulsive) potential. 

We believe that our approach represents the minimal 
model that allows one to study the dependence of the re¬ 
laxation rate on the potential strength and density with¬ 
out recourse to perturbative methods. 


III. NUMERICAL TECHNIQUE 


A. Solving the Schrodinger Equation 


For a numerical determination of the spin autocorre¬ 
lation function we solve explicitly the one-dimensional 
Schrodinger equation 


d t | t(>{x, s z ,t )) = -iH \tf)(x, s z , t )), (15) 


where ^(x, s z , t) is a nucleon single-particle wavefunction 
which depends on the one-dimensional coordinate a:, the 
z-conrponent of the spin s z = ±1/2, and time t. The 
Hamiltonian is 


H = — 


1 


_ 

2toat dx 2 


■Hi 


int • 


(16) 


This partial differential equation is solved by discretiz¬ 
ing it on a spatial grid of spacing a so that Xj = j a, 
j = 1, • • •, n x , i.e. the box size is L = ( n x + 1) a. We use 
Dirichlet boundary conditions ip( 0) = ip(L) = 0, imply¬ 
ing that there are n x grid points where ip{xj) ^ 0 and 
thus n x orthonormal states for each s z = ±1. For the 
kinetic energy we use the discretized prescription 


rzi „/.w_ ^ + V’fe-i) - 2 il>(xj) 

{-H-km W)\ x j) — 0 2 

2 nriNd 


(17) 


We then integrate over a time interval [0, t m ax] in steps 
of length dt. 

In practice, at each time step the complex linear equa¬ 
tion 

|^(t ± dt)) = \ip{t)) + Q + ^2 “ X ) ^ + dt ^ 

(18) 


is implicitly solved for \ip(t + dt)). This finite-difference 
representation of exp [—iH dt] is accurate to second order 
in dt and unitary. Equation (|l8|) is employed alternat- 
ingly for the n x -dimensional spatial part using 7/kin for 


both spin polarizations, and using the two-dimensional 
spin-dependent part Hi nt at all spatial grid points. 

In terms of a complete set of orthonormal functions 
| j, s z ) the correlation function is given by the expression 



(j, | e- H / T+iHt+ikx cj z e- lHt - ikx a z \j, s z ) 

3 Z 


where the partition function is 


(19) 


Z = 


j —1,... ,n x 

s z =± 1/2 


O'. 


,-H/T 


1 3,s z ). 


( 20 ) 


The outer average in Eq. © is over many configurations 
of the scattering potentials. 

The scalar products in Eq. (|19|) can be written 


as 


the matrix elements A,t{j, s z \e lkx a z \j, s z )B,t between the 


states | j,s z ) A ,t 


= e - iHt °- 


e H/T \ j,s z ) and \j,s z ) B ,t = 


e lHt e lkx a z \j,s z ). Put another way, for a complete set 


of basis functions | j, s z ) we need to solve the Schrodinger 
equation for the states e~ H ^\ j,s z ) and e~ lkx a z \j, s z ), 
respectively, and then take the matrix element of the op¬ 
erator e lkx a z . 

The runs presented here will always be for the long- 
wavelength limit with k = 0. We have verified that 
S(u>,k) is indeed nearly independent of k as long as 
k < T. Furthermore, we will always calculate the sym¬ 
metric (real) correlator R(t) = R eR(t). 


B. Long-Time Behavior 

The numerical solutions for R(t) obtained in this way 
have the somewhat unintuitive property that they relax 
to a nonvanishing value lim^oo R(t) = R oo > 0. This 
long-time behavior is an artifact of the discrete nature of 
our solution. 

This is understood most easily when we write the struc¬ 
ture function in terms of the energy eigenstates | n), en¬ 
ergy Wn, as 


S(ujk) = — “" /T K n l' 5 '(°’ k )l m >l 


n N 


3 Z 

x ± Lu n cj m ), (21) 


where Z = X^ n ( n l n ) e~ ul ' n ^ T is the partition function. 

For a finite number of states the weight of the diago¬ 
nal part of the double sum (ui n = to m ) does not vanish 
relative to the nondiagonal part (u> n ^ w m ). Therefore, 
5 , (w, k) inevitably has a nonvanishing amount of power 
proportional to S(u>) even though its spectrum otherwise 
is continuous because the ensemble average over different 
potential configurations washes out the discrete nature of 
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the energy differences to n ~ oj m for n fy m. The power of 
S(uj) shows up as a positive constant term in the Fourier 
transform R(t). 


C. Short-Time Behavior 

The relaxation behavior of R(t) cannot be exponential 
at arbitrarily early times, corresponding to arbitrarily 
large frequencies in Fourier space. Very high frequencies 
correspond to the (spatial) Fourier components of the 
scattering potential which is not arbitrarily hard. One 
expects that the (quantum) structure function has finite 
moments of the form 

, 22 ) 

for all n. (The case n = 1 is the usual /-sum.) They 
will exist if the potential is infinitely often differentiable 
everywhere. Assuming that this is the case we conclude 
that S(ui) must fall off faster than any power of u and 
that the correlation function can be represented as 

R(t) = 1 - (r 2 t ) 2 + e>[(r 2 t) 4 ], (23) 

i.e. it has a vanishing derivative at t = 0. 


D. Relaxation Rate 


The interpretation of the numerical R(t) would be 
straightforward if it were exactly an exponential decay 
law, which in general it is not. Therefore, we need a 
useful prescription to quantify the “relaxation rate,” or 
equivalently, the “width” of S(uj). Put another way, we 
need to define what we mean with T for a non-exponential 
behavior of R(t). 

Among many possibilities which differ by a factor of 
order unity we find it most natural, for the present pur¬ 
poses, to use the neutrino-nucleon cross-section reduc¬ 
tion as a starting point. By virtue of Eq. (B8) we can 
calculate the cross-section suppression from the numeri¬ 
cal R(t) : and then translate it into an equivalent width 
r eff of a Lorentzian that would produce the same cross- 
section reduction. From Fig. |l5] we conclude that all of 
the relevant information of R(t) is contained in the time 
interval 0 < t < 2/T. 


E. Born Approximation 


In one dimension for one scattering potential of the 
form Eq. 0 ) the spin-density structure function in Born 
approximation in the long-wavelength limit for to > 0 is 


OBorn(WJ = - 




dp ft 


V(p+) + v(p_) 


{p 2 + 2mArw) 1//2 

) 

dp f p (24) 


where V{k) = f+^dxe lkx V(x ) is the Fourier trans¬ 
form of the potential V(x), p± = ( p 2 + 2mjvw) 1 / 2 ± p , 
f p = e ~P~/( 2m NT) j g thg nuc leon occupation number, and 
p s is the (linear) density of scattering centers. This struc¬ 
ture function is not resummed, i.e. it diverges for to —> 0, 
as expected for a bremsstrahlung result without includ¬ 
ing multiple scattering. 

For a Gaussian potential of the form Eq. ( fi~3| ) the Born 
approximation should be valid when 

< 25 » 


where p is a typical nucleon momentum. 

If we take the Gaussian to be very narrow, approaching 
a 5-function, then the Fourier transform V is a constant 
and we find 


^Born ( w ) 


4(2tt) 1 / 2 Psm^ 2 (V 0 b) 2 r°° e-“ 2 

3 Ti/2 W 2 J o au J u 2 +u/t - 

(26) 


This expression diverges at ui = 0 because it is a first- 
order un-resummed perturbative result. We can still cal¬ 
culate the neutrino cross-section reduction by virtue of 
Eq. (B2), which in turn can be translated into an equiva¬ 
lent width r° orn of a Lorentzian which produces the same 
cross-section reduction. We find 


F° 

Born _ 

T 


1/2 

Psmjf 

T 3/2 


(V 0 b) 2 


(27) 


with 


C = 


dx 


1 — (1 + x + x 2 /6) e 


(27r) 1 / 2 Aa; 2 


3-7T 2 — 16 
18 (27T) 1 / 2 A 


= 0.68752, 


/ du — F == 

10 v« 2 + X 


(28) 


where the constant A is defined in Eq. (p6|). 


For sufficiently weak potentials V(x) a perturbative 
calculation of the structure function can be performed 
in Born approximation which is a useful test of our nu¬ 
merical code. For the 3-dimensional case and a single 
scattering center such a comparison has been performed 
in detail in Ref. [llj. 


F. Collision Rate 

We finally estimate the rate by which a nucleon in our 
model encounters a scattering potential. If the nucle¬ 
ons are at temperature T, their average kinetic energy in 
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one dimension is ^ mN{v 2 ) = \T so that ( v 2 ) = T/mjv- 
Taking a Maxwellian velocity distribution tells us that 
(|a[) = (2 /tt) 1 / 2 (ti 2 ) 1 / 2 so that (|r>|) = (2T /ttitin ) 1 ^ 2 • If 
the average (linear) density of scatterers in our model is 
p s , then our average collision rate is r co ii = (|n|)/? s , or 


r co ii = / 2 \ 1/2 

T \irTrriN J ^ 


(29) 


Ultimately, the purpose of our numerical runs is to check 
if this estimate can be reached, or even surpassed, for a 
plausible configuration. 


IV. NUMERICAL RESULTS 


B. Born-Approximation Test 

In order to test our numerical code we compare the 
(un-resummed) dynamical structure function in Born ap¬ 
proximation of Eq. ( |24| ) with the one derived from a nu¬ 
merical run. The characteristic parameters for the run 
are given in Table |. For these parameters, the condi¬ 
tion Eq. ( |j ) would be Vo < 35. Further, TBorn/7 1 = 
7.67 x KT^and T co n/T = 0.094, again implying that a 
perturbative estimate should be reasonable. 


TABLE I. Parameters for Run Bom Test 


Uo 

b 

RC“ 

T 

ps 

a 

L 

^max 

dt 

5 

1 

yes 

30 

0.1 

0.7 

100 

81.92 

0.005 


“Repulsive Core in the potential. 


A. Parameters 


We now turn to several series of numerical runs where 
we have calculated R(t). A given run is characterized 
by the maximum potential strength Vo (in MeV) and its 
width b (in fm) according to Eq. (jl||) and by whether or 
not we include a spin-independent repulsive core (RC) 
to avoid bound states. Further parameters are the tem¬ 
perature T (MeV), the grid spacing a (fm), the density 
of scatterers p s (fm -1 ), the box size L (fm), the total 
number of scatterers in the box n s , the number of config¬ 
urations calculated n c so that the statistical significance 
is proportional to n s n c , and finally the integration time 
W and the time step dt (both in MeV -1 ). We usually 
take the number of time steps t max /dt to be a power of 
2 as required for a numerical cosine transform of R(t) to 
determine S(co). 

As already indicated, the energy unit is MeV, the time 
unit MeV -1 , and the length unit 1 fm = (197 MeV) -1 . 
Note also that the distance between nucleons scales as 
r4 /3 = 77.3 p\{ 3 MeV and n„ 1/3 = 2.55 p^ /3 fm. We 
will give all quantities in these units in the remainder of 
this paper. Most runs presented in the following were 
performed for 5=1 and T = 30, for which Eqs. (p4|), 
(B2), and (|B5|) yield 


= 3.07x 10 -4 p s U o 2 . 


(30) 


Likewise, the estimated collision rate from Eq. (g9|) is 
numerically 


Tcoll 

~T~ 


5.13 


Ps 

T 1 / 2 ' 


(31) 


These parameters allow one a first estimate of what to 
expect in a given run. 

The required CPU-time scales roughly linearly with 
the step size t m&x /dt and with the number of configu¬ 
rations n c , but quadratically with the number of grid 
points L/a. 


In Fig. |I] we compare the numerical S(ui), with the 
Born approximation >§Born(w) = <S , Born(w) (1 + e -a, / T )/2 
where SBorn(w) is from Eq. (|24|). In addition, we show 
the Born approximation of Eq. (|20|) for a S- 

function potential with the same Vo b. The two Born 
curves differ only for large to where the structure of the 
potential is resolved in a collision. The numerical struc¬ 
ture function and the Born approximation agree very 
well for intermediate frequencies as they should. This 
agreement is also reflected in the effective spin fluctua¬ 
tion rate, T e ff = 0.0199, and Tbo™ = 0.0230. Note that 
Tgorn = 0.0487 is a factor 2.12 larger because a (5-function 
potential produces more power at large frequencies. 



FIG. 1. Numerical result for S of the Born Test run (solid 
line), as compared to ^Born (long dashed line). The dotted 
line is Sjon,, corresponding to a (5-function potential. 
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C. Variation of the Box Size 


In Sec. Ill B we have remarked that the numerical R(t) 
relaxes to a constant value for t —> oo, an effect attributed 
to the finite box size of our numerical runs. In order 
to verify that this constant term indeed decreases with 
increasing box size, we have performed a series of runs 
Box Size with the characteristics summarized in Table ||; 
there was no repulsive diagonal potential. 


TABLE II. Common parameters for Series Box Size. 


To 

b 

RC 

T 

a 

Ps 

^max 

dt 

30 

1 

no 

30 

i 

0.1 

32.768 

0.002 



TABLE III. 

Results for Series Box Size. 


Run“ 

n s 

n c 

Rex, 

1/Too 

B0 

i 

4096 

0.336 

2.98 

Bl 

2 

2048 

0.214 

4.67 

B2 

4 

1024 

0.157 

6.38 

B3 

8 

512 

0.127 

7.87 

B4 

16 

252 

0.105 

9.49 

B5 

32 

74 

0.104 

9.63 

B6 

64 

18 

0.090 

11.12 


“Run B?i has a box size L = 10 x 2". 


We show the numerical R(t) for several runs in Fig. 
and summarize the results in Table III. We have calcu¬ 
lated the limiting value for R(t) as the average of R(t) for 
20 < t < fmax- It becomes indeed smaller for an increas¬ 
ing box size, but only very slowly. Run BO with only one 
scattering center does not relax to a constant at all; the 
remaining oscillations of R(t) signal that “bound-state 
effects” in the small box are important. 



FIG. 2. Spin autocorrelation function R(t) for the Box Size 
runs BO, Bl, B2, B3, B6 (top to bottom). 


In Fig. U we plot 1/Roo against log 2 (n s ), revealing an 
almost linear behavior. The runs with the largest box 
size evidently still have large errors so that their exact 
behavior remains numerically uncertain. 

For the neutrino cross-section suppression, only the 
range t < 2/T « 0.07 is important (see Appendix B), 
where all R(t), even the one from run BO, agree perfectly 
with each other. Evidently, finite-box effects are of no 
relevance if the turn-over to the long-time plateau occurs 
at sufficiently late times, i.e. at t > 2/T. 

A similar conclusion is reached by studying the struc¬ 
ture function S(oj) which differs between different box 
sizes only for the few lowest Fourier modes (Fig. ||). Put 
another way, the power which goes into S'(O) is missing 
from the few lowest ui ^ 0 modes, but not from the high 
modes which remain unaffected by finite-box effects. 



FIG. 3. Asymptotic value Too = limt_,oo R{t) for the Box 
Size runs. 



FIG. 4. Lowest Fourier modes of R(t) for some of the Box 
Size runs. 
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D. Variation of the Potential Strength 


1 


When the scattering centers are “dilute” so that the 
nucleons can be thought of as travelling undisturbed be¬ 
tween collisions, the spin relaxation rate should satu¬ 
rate with increasing potential strength—the spin can be 
flipped at most once per encounter. To demonstrate this 
effect and to determine the saturation value of the po¬ 
tential we have performed a series of runs with a Gaus¬ 
sian potential of varying strength Vo. We have included 
a spin-independent repulsive potential to avoid bound 
states which would inevitably appear as the potential is 
turned up. The common parameters for this series Po¬ 
tential Strength are given in Table IV while parameters 
and results for specific runs are summarized in Table |v|. 
The runs are called P n with n the potential strength Vo 
in MeV. 



FIG. 5. Correlation function R(t) for some of the runs of 
the Potential Strength series. The curves are marked with Vq. 


TABLE IV. Parameters for series Potential Strength. 


b 

RC 

T 

a 

Ps 

n s 

dt 

r coii 

1 

yes 

30 

i 

0.1 

20 

0.002 

2.81 


TABLE V. 

Results for series 

Potential Strength. 

Run“ 

^max 

n c 

r Born 

Ti 

Roo 

Teff 

P2 

65.536 

22 

0.00368 

— 

— 

0.0033 

P4 

32.768 

31 

0.0147 

— 

— 

0.0135 

P8 

16.384 

135 

0.0589 

0.34 

0.10 

0.048 

P12 

16.384 

89 

0.132 

0.55 

0.15 

0.098 

P16 

16.384 

40 

0.235 

0.75 

0.17 

0.158 

P20 

8.192 

107 

0.368 

0.93 

0.20 

0.22 

P40 

8.192 

90 

1.47 

1.23 

0.27 

0.57 

P80 

8.192 

90 

5.89 

1.13 

0.33 

0.95 

P160 

8.192 

66 

23.5 

0.91 

0.37 

0.96 


“Run Pn has a potential strength Vo = n. 



FIG. 6. Short-time behavior of R(t) for the runs P20 and 
P80 of the Potential Strength series. 


In Fig. H we show R(t) for most of the runs; the curve 
for P16 falls almost exactly on top of P20. The runs 
P16-P160 show essentially the same relaxation behavior, 
except that the long-time value Rea increases with in¬ 
creasing potential strength. This is understood because 
for large potentials there are fewer low-lying energy levels 
that can be occupied by thermal nucleons, decreasing the 
effective number of energy eigenstates and thus increas¬ 
ing the importance of the diagonal term in the double 
sum of Eq. ©■ 

In order to quantify the overall relaxation behavior of 
the R(t) curves we approximate them by the ansatz 

R(t) = R 00 + (1-R 00 )e~ r ' t . (32) 

The best-fit parameters are given in Table |v|. For the 
runs P2 and P4 a global fit was not meaningful because 
the overall integration time was too short. 


The overall relaxation rate saturates and the satura¬ 
tion value is crudely estimated by F co ii- Still, with in¬ 
creasing potential the short-time behavior changes signif¬ 
icantly for the large potentials as illustrated in Fig. y, i.e. 
an increasing Vo still implies increasing power at higher 
modes of S(oj). However, R(t) falls off more slowly than 
a global exponential because of its vanishing derivative 
at t = 0. Hence the neutrino-cross section reduction is 
smaller than expected from the relaxation rate Ti, and 
saturates at larger values for Vo- The entries in Table |v| 
for r e ff confirm this trend; the saturation value for r e fi is 
very similar to that of Ti, i.e. about 1 MeV, but the po¬ 
tential strength at which this value is reached is Vo ~ 16 
for Ti, as compared to Vo — 80 for r e g. 

We may also compare r e fj with the simple Born es¬ 
timate. Evidently, the Born approximation gives us a 
reasonable estimate of r e fj and thus of the cross-section 
suppression up to Vo < 20, but substantially overesti¬ 
mates r eff for higher potential strengths. 
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Finally, we show in Fig. ^ the dynamical structure 
function S(ui) obtained from cosine-transforming R(t). 
We actually show u> 2 S(u), a quantity that would asymp¬ 
totically reach 2r for large u> if S(u) were a Lorentzian. 
For Vo > 12 the power spectra develop a conspicuous 
resonance feature. 



FIG. 7. Dynamical structure function for the runs of the 
Potential Strength series. 


Such resonance effects are even more important if 
we use a potential without a spin-independent repulsive 
term. This is illustrated by a similar series of runs with 
common characteristics sum mariz ed in Table VI and run- 
specific information in Table VII. Of course, in the per¬ 
turbative regime a spin-independent potential has no ef¬ 
fect on the spin relaxation rate so that the presence or 
absence of this term begins to show up only when the 
impact of the potentials is no longer small. The impact 
of the repulsive core on R(t) is illustrated by Fig. 0 for 
Vo = 20 and 160. The difference in T e ff is still ~ 10% at 
Vo = 5. 



FIG. 8. Correlation functions with (dotted) and without 
(solid) a repulsive spin-independent core for the potential. 


TABLE VI. Parameters for series Resonances. 


b 

RC 

T 

a 

Ps 

r coii 

1 

no 

30 

i 

0.1 

2.81 


TABLE VII. Results for series Resonances. 

Run“ 

^max 

10 3 dt 

n s 

n c 

r Born 

Left 

R5 

65.536 

2 

20 

23 

0.0230 

0.023 

R10 

32.768 

2 

20 

54 

0.0920 

0.091 

R20 

16.384 

2 

20 

92 

0.368 

0.37 

R30 

16.384 

2 

10 

240 

0.828 

0.83 

R40 

8.192 

1 

5 

932 

1.47 

1.51 

R60 

8.192 

0.5 

5 

265 

3.31 

3.5 

R80 

4.096 

0.25 

5 

280 

5.89 

6.3 

R120 

4.096 

0.25 

5 

301 

13.2 

14.8 

R160 

2.048 

0.125 

5 

213 

23.5 

28.5 


“Run Rn has a potential strength Vq = n. 


In Fig. we show the dynamical structure function 
for some of these runs. When bound states are im¬ 
portant, they have far more power at large frequencies 
than the previous runs with a repulsive core. Likewise, 
the effective relaxation rate r e ff increases with increasing 
Vo; there is no saturation. Indeed, both rBorn and r e ff 
scale as V 0 2 . For all runs rBorn/r e ff ~ 1, and thus the 
Born approximation is surprisingly good, despite being 
not expected to be applicable according to Eq. (^5|) . For 
Vo = 160 we reach r e ff/T = 1, but of course, such large 
potentials are far beyond what appears plausible for a su¬ 
pernova core where bound-state effects probably do not 
dominate. 



FIG. 9. Dynamical structure function for some of the runs 
of the Resonance series. 
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E. Variation of the Potential Width 


TABLE X. Parameters for series Density. 


In order to check to which degree the saturation value 
of r eff for large potential strengths Vo depends on the 
width b of the potential, we have calculated for a 
series of combination of Vo and b. For these runs the 
integration time was taken very short, just long enough 
to determine r e ff, but not long enough to calculate a de¬ 
tailed power-spectrum or other global information. The 
common param eters for this series Potential Width are 
given in Table VIII while the F e fj-values are shown in 


Table IX 


TABLE VIII. Parameters for series Potential Width. 


RC 

T 

a 

Ps 

n s 

L 

dt 

^max 

n c 

r coii 

yes 

30 

0.5 

0.1 

5 

50 

0.001 

0.1 

500 

2.81 



TABLE IX. 

r e ff for 

series 

Potential Width. 



25 

50 

75 

Vo 

100 

125 

150 

175 

200 

b = 0.5 

0.13 

0.33 

0.47 

0.55 

0.57 

0.58 

0.54 

0.52 

1.0 

0.31 

0.73 

0.96 

1.08 

1.07 

1.03 

0.95 

0.93 

1.5 

0.50 

1.12 

1.48 

1.63 

1.59 

1.54 

1.46 

1.39 

2.0 

0.68 

1.53 

2.03 

2.19 

2.12 

2.11 

1.98 

1.90 


The saturation value of T e ff increases with b , but for 
widths which represent the range of the nuclear force, it 
never reaches the “collision limit” T co n. 


Vo 

b 

T 

a 

L 

^max 

dt 

20 

1 

30 

0.5 

50 

4.096 

0.001 



TABLE XI. 

Results for Density (repulsive core). 


Run 

Ps 

n c 

r coii 

r Born 

r e fi 

D10 

0.2 

166 

5.6 

0.736 

0.46 

D20 

0.4 

164 

11.2 

1.47 

0.90 

D40 

0.8 

164 

22.5 

2.94 

1.76 

D60 

1.2 

164 

33.7 

4.72 

2.53 

D100 

2.0 

146 

56.2 

7.36 

4.1 

D200 

4.0 

145 

112 

14.7 

8.1 

D600 

12.0 

340 

337 

44.2 

27.9 


TABLE XII. 

Results for Density (no 

repulsive 

core). 

Run 

Ps 

n c 

r coii 

r Born 

Left 

El“ 

0.02 

56 

0.56 

0.0736 

0.067 

E3 6 

0.06 

60 

1.69 

0.221 

0.234 

E5 6 

0.10 

54 

2.81 

0.368 

0.38 

E7 6 

0.14 

60 

3.93 

0.515 

0.51 

E10 

0.2 

166 

5.6 

0.736 

0.72 

E20 

0.4 

166 

11.2 

1.47 

1.39 

E40 

0.8 

166 

22.5 

2.94 

2.69 

E60 

1.2 

166 

33.7 

4.42 

4.0 

E100 

2.0 

154 

56.2 

7.36 

6.4 

E200 

4.0 

163 

112 

14.7 

11.6 

E600 

12.0 

250 

337 

44.2 

40.0 

^max 

= 4 x 4.096 and dt = 

0.002. 




b t m ax = 2 x 4.096. 


F. Variation of the Density of Scatterers 


As a next test we have varied the density of scatterers, 
keeping the potential strength fixed at Vo = 20 where 
the efficiency of spin flips saturates according to the re¬ 
sults of the previous section. We have performed two 
series of runs, one (D) where the potentials include a re¬ 
pulsive core to avoid bound states, and one (E) without 
this term. The common parameters for both series are 
given in Table |x|, while parameters and results specific 
to the series D and E are given in Tables [Xl] and XII, 
respectively, and the dynamical spin structure functions 
for some of these cases are shown in Figs. and |lT 


The runs are named Dn or En where n is the number of 
scatterers in a given configuration. 

The effective relaxation rate grows linearly with den¬ 
sity (Fig. [if) even for densities so large that the poten¬ 
tials significantly overlap with each other. We have ob¬ 
served that this holds for any potential strength as long 
as most of the energy eigenfunctions are not localized. 
The absence of saturation can also be seen from the dy¬ 


namical structure functions themselves, Figs. [1(1 and 11 


TABLE XIII. Results for Density (repulsive core, Vo = 80). 


Run 

ps 

n c 

r coii 

r Born 

r e fr 

DV3“ 

0.06 

55 

1.69 

3.53 

0.609 

DV10 ft 

0.2 

40 

5.6 

11.8 

1.95 

DV40 C 

0.8 

47 

22.5 

47.1 

7.89 

DV200 C 

4.0 

425 

112 

235 

54.3 

DV600 C 

12.0 

108 

337 

706 

104 


“fmax = 2 x 4.096. 
b dt = 0.0005. 

Tmax = 4.096/2 and dt = 0.0005. 


Note that at high densities the structure functions ex¬ 
hibit the flat behavior at low frequencies expected from 
the classical picture, Eq. (jf), in the case without a re¬ 
pulsive core, whereas in the presence of a repulsive core 
they increase with frequency. 

As soon as the potentials are strong enough and/or 
the density is high enough for a significant fraction of the 
energy eigenfunctions to become localized we observe a 
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cj/MeV 

FIG. 10. S(u>) (upper panel) and u) 2 S(uj) (lower panel) for 

some of the runs of the Density series with repulsive core. In 
decreasing order at the left end, the lines represent the runs 
DIO, D40, D200, and D600. 



cj/MeV 

FIG. 11. Same as Fig. 0 but for the Density series with¬ 
out repulsive core., i.e., the runs E10, E40, E200, and E600. 


saturation in p Sl as demonstrated in Table XIII which 
lists results for potentials of strength Vq = 80 and with a 
repulsive core. This can be understood qualitatively by 
observing that |(n|<r(0, k)|n)| 2 is independent of the box 
size L for a localized eigenstate |n), as opposed to de¬ 
creasing with L for a non-localized state, and thus gives 
rise to an asymptotic non-vanishing Roo in the contin¬ 
uum limit L —> oo, preventing the scattering cross sec¬ 
tion (a) to decrease further [see Eq. (jB8|)]. At the same 
time, we expect the test-particle approximation adopted 
in the present work to break down in the regime with a 
significant fraction of localized states. 

Outside this regime, the highest effective spin fluctua¬ 
tion rate and the strongest scattering cross section sup¬ 
pression we found were T e g/T ~ 2 and (a)/<Jo — 0-5 for 



Vo — 80, Ps — 4, about a factor 2 smaller than r co u, see 
Tabic [Xlll. 

Finally, to get a better overview of the saturation be¬ 
havior for large potential strengths Vq in the presence 
of a repulsive core, we have calculated T e g for a series 
of combination of Vo and p s . Again, for these runs the 
integration time was taken very short. The common pa¬ 
rameters for this series Strength-Density are given hi Ta¬ 
ble XIV while the T e ff-values are shown in Table XV. 


TABLE XIV. Parameters for series Strength-Density. 


b 

RC 

T 

a 

L 

dt 

^max 

n c 

1 

yes 

30 

i 

50 

0.001 

0.1 

500 


TABLE XV. F e ff for series Strength-Density. 



25 

50 

Vo 

75 

100 

150 

ps = 0.1 

Fcoii = 2.81 

0.32 

0.70 

0.92 

1.03 

0.99 

0.5 

14.0 

1.59 

3.6 

4.7 

5.4 

5.2 

1.0 

28.1 

3.1 

7.1 

9.7 

10.9 

11.4 

2.0 

56.2 

6.2 

14.5 

22. 

27. 

38. 

4.0 

112.4 

12.8 

33. 

50. 

64. 

87. 


Again, even for extreme combinations of density and 
potential strength, does T e ff stay safely below the colli- 
sional estimate r co n. For combinations which may plau¬ 
sibly mimic nuclear matter, r e g always stays below T. 
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V. DISCUSSION AND SUMMARY 

We have numerically studied the relaxation behavior of 
a single nucleon which moves in a one-dimensional box, 
filled with a certain density of spin-dependent scatter¬ 
ed. This simple model is meant to mimic the collisional 
broadening of the dynamical nucleon spin-density struc¬ 
ture function in a supernova core. 

In most runs, the potential was modelled as a Gaussian 
with a width of 1 fm as suggested by the range of the 
one-pion exchange force which is the primary cause for 
nucleon spin fluctuations. However, the exact range of 
the potential has no significant impact on our conclusion. 

In order to quantify the spin relaxation behavior by a 
single number we have chosen an effective spin relaxation 
rate T e ff which is defined as the width of a Lorentzian 
structure function which produces the same neutrino- 
nucleon cross section suppression. For sufficiently weak 
potentials, our numerical runs reproduce the perturba¬ 
tive value in Born-approximation. 

When the potential has a repulsive core such that 
bound states cannot form, the spin-relaxation rate satu¬ 
rates with increasing potential strength at a value which 
is a factor 2-3 smaller than one would have expected 
from a classical analogy where one pictures the spin as 
being randomized each time the nucleon encounters a 
scattering center. The saturation value for the potential 
strength roughly coincides with the value when the Born 
approximation begins to break down, and where bound 
states begin to form in the absence of the repulsive core. 
For the range 1 fm of our Gaussian potential this is at a 
potential strength of around 20 MeV which also is typical 
for the spin-interaction energy between two nucleons. 

Furthermore, we observed that the spin-relaxation rate 
is roughly linear in the density of scatterers. This is ex¬ 
pected in the classical picture, but our simulations show 
that it even holds if the potentials overlap substantially. 
Saturation is only observed for densities high and/or po¬ 
tentials strong enough for the formation of a substantial 
fraction of localized energy eigenstates. However, in this 
regime, a strong correlation of neighboring nucleon spins 
is expected, and the test-particle approximation adopted 
in our work is no longer meaningful. 

When the scattering potentials lead to bound states 
(no repulsive core), there is no saturation of r e ff, not with 
increasing potential strength and not with an increasing 
density of scatterers. Instead, the spin-relaxation rate 
follows the Born approximation amazingly closely, de¬ 
spite the fact that it is not applicable. However, the 
test-particle approximation again is not viable as a model 
for a nuclear medium because the bound states which 
dominate in this regime signal the importance of spin- 
spin-correlations. 

In any case, we expect the case of repulsive potentials 
without strongly bound states to be more appropriate as 
a model for the conditions in a supernova core because of 
the relevant properties of nuclear matter. In this situa¬ 


tion, the value r co n derived from a classical picture pro¬ 
vides an upper limit for the spin-fluctuation rate which 
saturates for potential strengths expected from nuclear 
interactions within about a factor 2, and captures the 
linear dependence on the one-dimensional density in the 
regime where spatial spin correlations are small and the 
test-particle approximation is applicable. 

For combinations of density and potential strength 
which do not look like absurd overestimates of what is 
appropriate for nuclear matter, the spin-fluctuation rate 
is always safely below the temperature, and the scatter¬ 
ing cross-section suppression is at most a few ten percent. 
The Born approximation, in contrast, significantly over¬ 
estimates the spin-fluctuation rate and predicts far too 
large cross-section reductions. 

Naturally, the conclusions from our one-dimensional 
model do not directly carry over to three dimensions, but 
can only be taken as an indication of what one might ex¬ 
pect. If the classical collision rate for three dimensions 
embodied by our Eq. (^o|) is taken as an upper limit 
for the spin relaxation rate, then again a naive Born- 
approximation calculation is a significant overestimate, 
and the true collisional broadening of S(ui, k) is too small 
to cause a significant cross-section reduction. 

Our calculations also indicate that the Born approxi¬ 
mation, while being a significant overestimate of the spin- 
relaxation rate, it is not an absurd overestimate. For the 
relevant conditions it may overestimate F by a factor of 
a few, but not by orders of magnitude. Therefore, we do 
not expect that quantities like the axion or neutrino-pair 
emission rate, which depend on the collisional broaden¬ 
ing of S(uj,k), suffer much of an additional suppression 
beyond those discussed in Ref. j|]. 

From our simple toy model we find no evidence for 
extreme deviations from the overall picture developed in 
Ref. (Hu]], which may be summarized by saying that 
spin fluctuations look small enough to avoid excessive 
neutrino opacity reductions, but large enough to avoid 
huge reductions of the axion or neutrino pair emissivities. 

Clearly, more accurate and detailed calculations should 
be performed in three dimensions with realistic nucleon 
interaction potentials. Furthermore, a treatment of the 
high-density limit should go beyond the test-particle ap¬ 
proximation and apply a full quantum Monte Carlo ap¬ 
proach @ to finite temperature and nuclear matter. 
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APPENDIX A: 

PERTURBATIVE ESTIMATE OF T 



V^al n N 
T !/ 2 m ^ 2 


1.24 puT 30 1/2 , 


(A3) 


From a perturbative calculation of the bremsstrahlung 
process NN —> NNvV one can easily extract the (unre- 
summed) dynamical structure function. For nondegen¬ 
erate nucleons and on the basis of a one-pion exchange 
potential it has been evaluated in Ref. 0; even for the 
conditions of a supernova core this potential should be 
a reasonable representation of the tensor force | 8 |. The 
result was written in the form 

S«M = ^s(aVT), (Al) 

where s(x) is a slowly varying dimensionless function of 
order 1, normalized such that s(0) = 1. Further, = 
4^/n a 2 nNT 1/,2 m N 5 ^ 2 with un the nucleon density, ton 
the nucleon mass, and a n = (/2 to/to w ) 2 /47 t ft 15 with 
/ft 1 the pion-nucleon “fine-structure constant.” The 
function s is explicitly 


where pu = p/10 14 g cm -3 and T 30 = T /30 MeV. Even 
if 7 is only a few is enough to imply a significant cross- 
section suppression (Appendix B). 


APPENDIX B: CROSS-SECTION REDUCTION 


The reduction of the axial-current neutrino-nucleon 
scattering cross section a by spin fluctuations was dis¬ 
cussed in detail in Ref. |7j]. With a 0 the vacuum cross 
section and (• • •) an average over a Maxwell-Boltzmann 
distribution of neutrino energies, one finds |7j 




du> 

2tt 


S(w) 


UJ 

T 




e"“ /T . (Bl) 


Equivalently, the amount of reduction can be expressed 
in the form 


s{x,y) 


e~ x / 2 + e x / 2 f°° t/2 3a; 2 + 6 ty + 5 y 2 

16 i\ x \ 1 6 3 

2 \Jt 2 — x 2 

^ _ 

x 2 + 2 ty + y 2 


1 ] no . ( * + V + Vt 2 - a; 2 ^ 
t + y gyt + y — \/t 2 — x 2 J 


(A2) 


where x = ui/T and pion-mass effects are included by 
virtue of the parameter y = m 2 /rriNT = 19.4 MeV/T. 
We show s(x, y) in Fig. |l^ for several values of y. 



FIG. 13. Dimensionless structure function according to 
Eq. (A2) for y = 0,..., 2 in steps of 0.5. 


For temperatures between 10 and 40 MeV, correspond¬ 
ing to y = 0.5-2, this function is approximately constant 
for a large range of x = ui/T. If we approximate it by 
1/2 we may write S^\u>) ft 2r/o ; 2 with T = r CT /4 or 


= (B2) 

where the dimensionless phase-space weight function is 
G{x) = 1 — (1 + x + x 2 /6) e~ x . (B3) 


This representation is needed to determine the suppres¬ 
sion when the structure function S(u) is given by an 
un-resummed perturbative expression and thus diverges 
at qj = 0 . 

We show the cross-section suppression Eq. (Bl) in 
Fig.y as a function of 7 = T/T for our generic structure 
function 


S(w) = 


2 r 


+ r 2 1 +, 


o/T 


implied by Eqs. ( 0 ) and (^j). Limiting cases are 

(< 7 ) _ f 1 - A'f for 7 < 1 , 

( 00 ) “ t B h for 7 > 1. 


(B4) 


(B5) 


The coefficients are 


A = 


dx 


e x - (1 + a; + a; 2 / 6 ) 


7T JO 


X 2 (l 


= 0.438703, 


B = — [ dx 2 + x + x2 / 6 = 1.59745. (B 6 ) 

7T 7 0 1 + e x 

We will need to know the value of 7 corresponding to 
a given cross-section suppression s = |<5(cr)/(cr 0 ) |. An 
approximation formula, good to within about 1 %, is 


7 («) = 


B (1 — y/l — S ) , 
1 — s 


AV 


(B7) 


It is essentially an interpolation between the inversions 
of the limiting functions of Eq. (B5). 
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r = r /t 

FIG. 14. Suppression of the average axial-current neu¬ 
trino-nucleon scattering cross section according to Eq. (Bl). 


In our present study, the primary quantity produced by 
the num erical runs is R(t). Therefore, we finally express 
Eq. (Bl) in its Fourier transformed version, 


(°) 

Wo) 


dr R(t/T ) F(t) 


(B8) 
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where 

^ 2 f°° J , , 2 + x + x 2 /6 

F(t) = - dx cos(zr)-——-, (B9) 

7r Jo i + e 

shown in Fig. [l^. It is normalized as / 0 °° dr F{t) = 1 
and F(0) =B = [tt 2 + 121og(4) + 3C(3)]/6 tt = 1.59745. 
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